Magic Number 7 ± 2 in Globally Coupled Dynamical Systems? 
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The prevalence of Milnor attractors has recently been reported in a class of high-dimensional 
dynamical systems. We study how this prevalence depends on the number of degrees of freedom 
by using a globally coupled map and show that the basin fraction of Milnor attractors increases 
drastically around 5-10 degrees of freedom, saturating for higher numbers of degrees of freedom. It 
is argued that this dominance of Milnor attractors in the basin arises from a combinatorial explosion 
of the basin boundaries. In addition, the dominance is also found in a system without permutation 
symmetry, i,e., a coupled dynamical system of non-identical elements. Possible relevance to the 
magic number 7 ± 2 in psychology is discussed. 

PACS numbers: 05.45 Jn, 87.19 La, 05.45 Ra 



In recent studies of dynamical systems with many de- 
grees of freedom, the prevalence of Milnor attractors has 
been recognized |l[ g]. The Milnor attractors are defined 
as follows: an arbitrary small perturbation to an orbit at 
a Milnor attractor can kick the orbit away from it to a 
different attractor, even though a finite measure of ini- 
tial conditions is attracted to the attractor by temporal 
evolution||, |, | . In other words, the basin of the attrac- 
tor touches the attractor itself somewhere. An orbit is 
often attracted to a Milnor attractor, but can be kicked 
away from it by infinitesimal perturbation. 

It should be noted that Milnor attractors can ex- 
ist in low dimensional dynamical systems like a two- 
dimensional map as well. When changing the parameter 
of a dynamical system, the basin boundary of an attrac- 
tor may move until, for a specific value of the parameter, 
the basin boundary touches the attractor. If, for this 
parameter value, the attractor has a positive measure of 
initial conditions forming the basin of attraction, it be- 
comes a Milnor attractor. It is naively expected that the 
above situation occurs only for very specific parameter 
values (e.g., at a bifurcation point), and that the Milnor 
attractors may not exist for an open set of parameter 
values. 

In recent studies of a class of dynamical systems, how- 
ever, Milnor attractors are found to be sometime preva- 
lent, occurring not only for specific isolated parameter 
values, but for continuous ranges of parameter values. 
Furthermore, the measure of the initial conditions that 
belong to the basin of these Milnor attractors is a rela- 
tively large proportion of the phase space. Indeed, for 
some parameter ranges, almost all randomly chosen ini- 
tial conditions fall onto Milnor attractors R 0. 

Such dominance of Milnor attractors is often found in 
high-dimensional dynamical systems, for example, glob- 
ally coupled maps with 10 degrees of freedom or so, and 
within a range of parameter values where many attractors 
coexist. The question we address in the present Letter is 
why can there be so many Milnor attractors in a "high- 



dimensional" dynamical system, and how many elements 
are sufficient for constituting such 'high' dimensionality. 
With the help of numerical results obtained from simu- 
lations of globally coupled maps, we will show that the 
dominance starts to be common at 5 ~ 10 degrees of free- 
dom. We propose a possible origin for this dominance of 
Milnor attractors, by noting a combinatorial explosion of 
saddle points (or more generally basin boundary points). 
We will also show that the prevalence of Milnor attrac- 
tors is observed even in a system without symmetry, i.e., 
in a coupled dynamical system of non-identical elements. 
Finally we briefly discuss the possible relevance of our re- 
sults to the magic number 7 ± 2 discussed in psychology. 

As a prototype example to study this problem we use 
a globally coupled mapB 



x n+1 (i) = (1 - e)f(x n (%)) + ^ £/(*»(?')). (1) 



where n is the discrete time and i the index for its ele- 
ments (i = 1, 2, • • - ,iV = dimension of the system). For 
the elements we choose f(x) = 1 — ax 2 , since the model 
has thoroughly been investigated as a prototype model 
for high-dimensional dynamical systems. The coupling 
parameter e is fixed at 0.1, since for this value the typ- 
ical behaviors of the above GCM that are relevant here 
can be observed by changing only a. 

In the present model, each attractor can be coded by 
the so-called clustering condition, that is to say, by the 
way how the N elements of the system partition into 
mutually synchronized clusters, i.e., a set of elements in 
which x„(i) — x n (j) B. Each attractor is coded by the 
number of clusters k and the number of elements Nk 
in each cluster with the clustering condition given by 
(Ni,N 2 , • • • , i\T fe )§. Indeed, for a GCM with N = 100, 
that many initial conditions are attracted to a Milnor 
attractor. Here, in order to discuss the dependence of 
the size of the basin fraction of the Milnor attractors on 
the number of degrees of freedom, we have computed the 
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FIG. 1: The basin fraction of Milnor attractors plotted as 
a function of the parameter a, for N = 3,5, 7, and 9. For 
the present simulations, we take 1000 randomly chosen initial 
conditions, and iterate 10 5 steps. Then the orbit is perturbed 
as x n (i) + 10~ <r, with a as a random number over [-.5, .5]. 
With 1000 trials of such perturbations, we checked whether 
the orbit remains on the same attractor or not. If some of the 
100 trials result in an escape from the original attractor, it is 
regarded as a Milnor attractor Q. 

ratio of initial conditions that are attracted to Milnor 
attractors. 

As shown in Fig.l, the basin fraction of Milnor attrac- 
tors is large around a w 1.65, especially for N > 5. Its 
parameter dependence, however, is quite strongly depen- 
dent on N, at least when N is not so large. Hence, it is 
not that relevant to compare the behaviors for different 
TV at a given value of a. Instead, we compute the average 
basin fraction of Milnor attractors over the parameter in- 
terval 1.55 < a < 1.72. In Fig. 2, this average is plotted 
as a function of the number of degrees of freedom N. The 
increase of the average basin fraction of Milnor attractors 
with ./V is clearly visible for N w (5 ~ 10), while it levels 
off for N > 10[|. 

Now we discuss how the dominance of Milnor attrac- 
tors appears. In a system with identical elements, due 
to the symmetry, there are at least M(N\, ■ ■ ■ , Nk) — 
H" m Ni \ noversctsofjv^jv, attractors for each cluster- 
ing condition, where mi is the number of clusters with 
the same value Nj |Q . Then, combinatorial increase in the 
number of attractors can be expected when many of the 
clustering conditions are allowed as attractors. However, 
the increase of the number of attractors cannot explain 
the increase of the basin for Milnor attractors. In Fig. 3. 
we have plotted the number of attractors and the basin 
fraction of Milnor attractors for N = 10. As can clearly 
be seen, the dominance of the Milnor attractors is not 
necessarily observed when the number of attractors is 
high. Rather, the fraction of Milnor attractors gets large 
even when many attractors start to disappear with the 
increase of a. 



FIG. 2: The average fraction of the basin ratio of Milnor 
attractors. After the basin fraction of Milnor attractor is 
computed as in Fig.l, the average of the ratios for parameter 
values a — 1.550, 1.552, 1.554, • • ■ 1.72 is taken. This average 
fraction is plotted as a function of N. (The dotted line depicts 
the case where a is incremented by 0.01 instead of 0.002; For 
N > 10, the simulation is carried out only with this choice.) 
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FIG. 3: The number of attractors (+) estimated from simu- 
lations over 10° initial conditions. The estimated number of 
attractors is plotted as a function of a. N — 10. All attrac- 
tors that are concluded to exist by the symmetry argument 
are also counted. The basin fraction of Milnor attractors ob- 
tained in the same way as in fig.l is also by dotted line with 
x. 

For the parameter region where many attractors start 
to disappear, there remain basin boundary points sepa- 
rating such (collapsed) attractors and the remaining at- 
tractors. To explain the prevalence of Milnor attractors, 
we discuss how the distance between an attractor and 
its basin boundary changes with N. Consider an inter- 
val having one basin boundary for x\ > x* and x\ < x* 
separating two attractors. If an iV-dimensional phase 
space is partitioned by (basin) boundaries at xi = x* 
for i = 1, ■ ■ ■ N, it is partitioned into 2 N units. In this 
case, the distance between each attractor and the basin 
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boundary does not change with N. 

On the other hand, consider a boundary given by 
some condition for (xi, • ■ • , Xn). In the present system 
with global (all-to-all) couplings, many of permutational 
change of Xi in the condition give also basin bound- 
aries. Here the condition for the basin is typically re- 
lated with some attractor (or a ruin of collapsed attrac- 
tor) that has clustering (N\, ■ ■ ■ , Nf.). Then there are 
M(Ni,---,Nk) partitions by boundaries equivalent by 
permutations. The number of regions partitioned by the 
boundaries increases combinatorially with N. Roughly 
speaking, it increases in the order of (N — 1)!, when the 
boundary has a variety of clusterings (i.e., large M). Re- 
calling that the distance between an attractor and the 
basin boundary remain at the same order for the parti- 
tion of the order of 2 N , the distance should decrease if 
(N — 1)! is larger than 2 . Since the former increases 
remarkably faster for N > 5 than the latter, the dis- 
tance should decrease drastically for N > 5. Since this 
argument is applied for any basin boundary condition 
with a complex clustering having combinatorially large 
M(Ni, ■ ■ ■ ,Nk), the distance between an attractor and 
such basin boundary may often go drastically small. For 
a certain parameter regime ( 1.64 < a < 1.68 in the 
present case), basin boundaries with such partitions are 
dominant, and the above argument is applied to most at- 
tractors. Although this explanation may be rather rough, 
it gives a hint to why Milnor attractors are so dominant 
when N is larger than 5 ~ 10. 

Since the above discussion is based mainly on simple 
combinatorial arguments, what we need is instability in 
orbits leading to many attractors, and global (all-to-all) 
coupling to allow for permutations of elements. Then, 

the dominance of Milnor attractors for N « (5 ~ 10) 
may be rather common in globally coupled dynamical 
systems. To check this possibility, we have also made 
numerical simulations for Josephson junction series ar- 
rays that are globally coupled through a resistive shunt- 
ing load and driven by an rf bias current pi, given by 



\-sin(f)j+ga/NY^ m=1 4>m = idc + i^[sin{Vl r ft). 
By using the parameter values adopted in Q, we have 
computed the basin volumes for Milnor attractors, at the 
partially ordered phase where a variety of attractors with 
many clusters coexists. Again, the basin volumes are 
close to for N < 4, and increase at 5 < N < 10, It is 
also interesting that pulse-coupled oscillators with global 
coupling also show the prevalence of Milnor attractors 
for N > 5 §. 

So far reports on the prevalence of Milnor attractors 
have been limited to a system with symmetry. For ex- 
ample, in the GCM (1), the permutation symmetry aris- 
ing from employing identical elements leads to a combi- 
natorial explosion in the number of attractors as men- 
tioned. Then, one may wonder whether the prevalence 
of Milnor attractors is possible only for such highly 
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FIG. 4: The basin fraction of Milnor attractors for a GCM 
with inhomogeneous parameters. Here a(i) = ao + .1 x (i — 
1)/{N — 1). Plotted as a function of the parameter ao, for 
N = 4 — 10. (For N = 3, the fraction is almost zero for all 
ao.) The fraction is computed in the same way as in Fig.l, 
except that 100 trials of perturbations were used instead of 
1000. Since the clustering condition cannot be used in this 
case, we checked whether the orbit is on the same attractor, 
by computing the temporal average of x n {i) over 5 x 10 6 steps, 
before and after each perturbation. If the average agrees to 
within a precision of 10 -3 , the orbit is regarded to be on the 
same attractor. (We have investigated other average times in 
the range from 10 4 to 10 5 with precisions ranging from 10 -4 
to 10~ 2 . Even though the basin fractions may change a little, 
the existence of the Milnor attractors and their dominance 
around ao = 1.61 for N = 7, 9 is invariant.) 



symmetric systems, especially because some Milnor at- 
tractors are known to disappear when introducing tiny 
asymmetries jllf . We have therefore studied a GCM 
with inhomogeneous parameters fl^| , given by x n +i(i) = 
(1 - e)fi(x n (i)) + J2j fj(x n (j))- with fi(x) = 1 - a,x 2 , 
and a, = a + a w j^. 

In Fig. 4, we have plotted the basin fraction for Milnor 
attractors with the change of the parameter ao while fix- 
ing a w = 0.1, for N = 3,4,- • • , 10. Although the fraction 
is smaller than in the homogeneous case, Milnor attrac- 
tors are again observed and found to dominate the basin 
volume for some parameter region. As in the symmet- 
ric case, the basin fraction of Milnor attractors increases 
around N « (5 ~ 10). 

Note that even though complete synchronization be- 
tween two elements is lost, clusterings as with regards 
to the phase relationships can exist |D3]. For example, 
there are two groups when considering the oscillations 
of phases as large-small-large... and small-large-small..., 
that are preserved in time for many attractors. Simi- 
larly, it is natural to expect an explosion in the number 
of the basin boundaries for some parameter regime. Ac- 
cordingly the argument on the dominance of Milnor at- 
tractors for a homogeneous GCM can be applied here to 
some degree as well. 
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The term magic number 7 ± 2 was originally coined in 
psychology when it was found that the number of 
chunks (items) that is memorized in short term memory 
is limited to 7 ± 2. (See [|l5] for a pioneering study to re- 
late this number with chaos.) To memorize k chunks of 
information including their order (e.g., a phone number 
of k digits) within a dynamical system, let us assign each 
memorized state to an attractor of a fc-dimensional dy- 
namical system, as is generally adopted in neural network 
studies. In this k dimensional phase space, a combina- 
torial variety of attractors has to be presumed in order 
to assure a sufficient variety of memories. Depending on 
the initial condition (given by inputs) , an orbit has to be 
separated to different attractors. Then, a combinatorial 
explosion of basin boundaries is generally expected with 
the increase of k, if the neural dynamics in concern are 
globally coupled (as often adopted in neural networks). 
Then, following the argument in the present paper, Mil- 
nor attractors may be dominant for k > (5 ~ 10). (Recall 
that the number does not strongly depend on the choice 
of models, since it is given by the combinatorial argu- 
ment.) Since the state represented by a Milnor attractor 
is kicked out by tiny perturbations, robust memory may 
not be possible for information that contains more than 
7 ± 2[jl6| chunks. Although this explanation is a rough 
sketch, it can possibly be applied to other systems that 
adopt attractors as memory [[il)). 

In dynamical systems, it is well known that the di- 
mensional cutoff > 3 plays an important role for the 
existence of chaos. It is interesting then to investigate 
whether there are certain higher dimensions that simi- 
larly form dimensional boundaries beyond which the be- 
havior of a dynamical system changes qualitatively. The 
present study may shed new light on this possibility. 
Also, it is interesting to note that in Hamiltonian dynam- 
ical systems, agreement with thermodynamic behavior is 
often observed only for degrees of freedom higher than 
5 ~ 10|2(J. Considering the combinatorial complexity 
woven by all the possible Arnold webs (that hence may 
be termed "Arnold spaghetti"), the entire phase space 
volume that expands only exponentially with the number 
of degrees of freedom may be covered by webs, resulting 
in uniformly chaotic behavior. If this argument holds, 
the degrees of freedom required for thermodynamic be- 
havior can also be discussed along the line of the present 
Letter. 
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